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Computation of Atomic Energy Levels: Spectrum of Singly- 
Ionized Tantalum (Ta n) 

Richard E. Trees, W. F. Cahill, and P. Rabinowitz l 

Energies and wave functions of atomic systems can be calculated as the characteristic 
values and vectors of matrices in accordance with long-established procedures. Coding for 
carrying out these computations on the Standards Eastern Automatic Computer is described. 
Some preliminary values are given for the energy levels of the spectrum of singly-ionized 
tantalum (Ta n). 



1. Introduction 

Calculations of atomic energy levels have been 
helpful in expediting the analysis of experimental 
data. 2 3 However, except in the simplest cases, the 
numerical work is excessive for hand computation, 
and the calculations have had very limited applica- 
bility. It is possible to overcome this limitation by 
utilizing digital computers. A calculation of the 
levels of the low even configurations in the spectrum 
of singly-ionized tantalum (Ta n) is given. This is 
a typical example of a calculation that is too difficult 
for hand computation, although fairly simple as 
compared with others in the field of complex spectra. 

2. Statement of Problem 

Spectra in which configuration interaction and 
spin-orbit coupling are both important are consid- 
ered; the basic theory for these spectra is well estab- 
lished. 4 5 According to this theory, energies of levels 
with the same total angular momentum (i. e., the 
same "J-value") are the eigenvalues of a matrix, 
which is here designated by {a tj ). The J-value is 
omitted for brevity, but there will be one of these 
matrices for each J- value that is allowed, in the con- 
figurations of the spectrum under consideration. 
The configurations occurring in Ta n are d*, d? s, 
and d 2 s 2 ; these give rise to matrices of orders 9X9, 
12X12, 21X21, 15X15, 15X15, 6X6, and 3X3 for 
the allowed J-values 0, 1, 2, 3, 4, 5, and 6, respec- 
tively. 

Each matrix element is a linear combination of a 
limited number of "parameters", p ik \ multiplied by 
known "coefficients" cf- (the parameters are the 
same for all J-values). 



a ii — 2-JkCii V 



(1) 



i Present address: Weizmann Institute of Science, Rehovot, Israel. 

2 R. E. Trees and M. M. Harvey, J. Research NBS 49, 397 (1952) RP2378 

3 M. Gehatiah, Phys. Rev. 94, 618 (1954). 

W. F. Meggers has recently located the c 2 D(3/2) level of Hf n at 30595, less 
than 100 K away from Gehatiah's calculated position of 30692. 

< E. U. Condon and G. H.Shortley, The theory of atomic spectra (Cambridge 
University Press, London, 1951). 

« G. Racah, Phys. Rev. 62, 438 (1942); 63, 367 (1943); 76, 1352 (1949). Referred 
to as R II, R III, and R IV, respectively. 



These linear forms can be derived in a straightfor- 
ward manner by use of methods developed by 
Racah (see footnote 5). For commonly occurring 
configurations, the calculation of any one element of 
the "coefficient matrix' ' (cJJ } ) involves at most the 
multiplication of a few factors, which are given in 
readily available tables. 6 The parameter associated 
with a given coefficient matrix is defined in theory as 
a radial integral (Slater integral, spin-orbit integral, 
etc.) (see footnote 4). This can be evaluated if 
radial wave functions are determined by the method 
of self-consistent fields, but the work involved in this 
procedure is great. In the present calculations, the 
parameters are adjusted by least-squares to get best 
agreement with levels that have already been identi- 
fied experimentally. This limits the use of the theory 
to spectra in which the analysis of the experimental 
data is already well started. This limitation may 
also be removed when computers are used to carry 
out self-consistent field calculations for complex 
spectra such as Ta n. 

The coefficient matrices for the electrostatic inter- 
action in the configurations of Ta n arc given in 
R ii and R in. For each J-value there are 16 of these 
associated with parameters identified by the letters 
A, B, C, G 2 , H 2 , and a (see table 2). 7 In addition to 
these, there are three parameters f associated with 
the coefficient matrices of the spin-orbit interaction. 
The latter coefficient matrices have been given for 
the configurations d 2 s 2 (see footnote 4) and for d* 
also. 8 The matrices for d 3 s have been computed by 
Trees; for the sake of brevity, the results are not 
given here. 

Recapitulating, there are 19 different parameters, 
and for each of the seven J-values there are 19 co- 
efficient matrices, one associated with each param- 
eter. 9 Though the orders of the coefficient matrices 

6 Because there are often many elements, it will nevertheless require much time 
and effort to set up these coefficient matrices correctly; probably automatic com- 
puters can be utilized in this phase of the work. 

7 For a discussion of the parameter « and other advances in atomic theory not 
covered in footnotes 4 and 5, see G. Racah, Kungl. Fysiografiska Sallskapets 
Handlingar N. F. 65, 31 (1955). (Proceedings of the Rydberg Centennial Con- 
ference on Atomic Spectroscopy). 

8 H. Greyber, Ph. D. thesis, University of Pennsylvania, 1953 (unpublished). 

9 It is convenient to use the code of section 6 to calculate the ^-values; to do this 
the Z/S-coupling ^-values are included as an extra coefficient matrix, the "par- 
ameter" associated with the <?-value always being zero. Therefore, 20 parameters 
were actually used for Ta n. 
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are the same as those already given for the appro- 
priate J-value, there are very few nonzero elements 
in most of the coefficient matrices; the only coefficient 
matrices with large numbers of nonzero elements 
are those associated with the two parameters f (rf 4 ) 
and f(d 3 s). 

3. Eigenvalues and Eigenvectors of Real 

Symmetric Matrices 

Preliminary estimates of the parameters p {k) are 
made by using methods given in section 7.1., and the 
matrix elements are then calculated according to 
formula (1). This can be done by hand computa- 
tion; or the coding described in section 5 can be 
used to generate the matrices on the National Bureau 
of Standards Eastern Automatic Computer (SEAC). 
The eigenvalues and eigenvectors of the matrices are 
then to be determined. 

A SEAC code that will do this is based on the 
Duncan-Collar iteration procedure. 10 The use of 
this code with matrices up to order 12X12 has been 
briefly reported. 11 The code can handle matrices 
with orders up to 25X25. The major disadvantage 
of this code is that it requires too much machine 
time for regular use when all eigenvalues of several 
high-order matrices are desired. 12 Because of this, 
no effort was made to determine all the eigenvectors 
and eigenvalues for the matrices of Ta n set up with 
the preliminary values of the parameters. No 
matrices were set up with the final values obtained 
for the parameters. 13 

4. Least-Squares Adjustment of Parameters 

The matrices set up with preliminary values of the 
parameters p (k) have the eigenvalues \ (m) and the 
eigenvectors X im) determined (the components of a 
vector will be indicated by x} m) ) as described in 
section 3. The eigenvalues are approximated as a 
function of the parameters by the following linear 
formula, which will be accurate to first order, in the 
difference between the parameters and their pre- 
liminary values: 



where 



X< w =S*d|5 ) P cw , 



J (m) % A >v ( m ) /> (k) /y. (m) 



(2) 
(3) 



io The matrix is modified to obtain successive values by the method given by 
W. M. Kincaid, Quart. Appl. Math. 5, 320 (1947). 

» D. E. Mann, L. Fano, W. F. Cahill, and T. Shimanouchi, J. Chem. Phys. 
33, 764 (1954). 

^Appreciation is expressed to W. G. Hall for running SEAC to get the energy 
values in Ta n reported in section 7.2. He also obtained values in the first spec- 
trum of tungsten (W i) for G. Racah; a report on the latter work, which involved 
matrices with orders up to 25 X 25, has been given in footnote 7. On the basis 
of his work, Hall has estimated that the average time required to iterate for one 
eigenvalue and eigenvector is about 5 min for a 10 X 10 matrix, 20 min for a 20 X 20 
matrix, and 50 min for a 25 X 25 matrix. 

is K. Goldberg has a code that determines only the eigenvalues; this code is 
similar to that described by W. Givens, ORNL 1574 (1954). It can determine 
the eigenvalues of matrices up to order 30 X 30; the time required for a 30 X 30 
matrix is 25 min. About 10 min are needed for a 20 X 20 matrix. This code is 
being extended to give the vectors; it is expected that the extended code will be 
able to handle matrices of orders up to 23 X 23. 



The coefficients c/*? are the ones already used in (1). 
The linear formulas (2) can be set up by hand; or 
the coding described in section 6 can be used to 
generate them on SEAC. A least-squares adjust- 
ment can then be made to find new values of the 
parameters p {k \ which will make the linear forms (2) 
agree best with a corresponding set of observed 
experimental values. 

A general-purpose orthonormalizing code has been 
set up for SEAC, 14 and it can be used to carry out 
least-squares calculations. The coefficients rfj^j, for 
fixed "&", are regarded as components of a vector in 
a space with the same dimensionality as the number 
of observations to be fitted. For successive ^-values, 
each vector is orthogonalized to all preceding vectors 
with smaller ^-values. The observed values, re- 
garded as an additional vector, are then expanded in 
terms of this set, and the residues that remain after 
the expansion will be the least-squares deviations. 
The process can be interrupted before the ortho- 
normalization is completed for all ^-values, so that 
the least-squares calculation can be carried out, say, 
first by omitting and then by including the last 
parameter, without making two separate calcula- 
tions. As many sets of observations may be fitted as 
desired in a single calculation. Consequently, differ- 
ent possibilities for correlating observations with 
theory may be tested. An observation may be 
omitted from the calculation by inserting a new 
"weight vector" in which the component correspond- 
ing to the deleted observation is given the weight 
zero. The vectors must be "augmented" in these 
applications (see footnote 14), so that the least- 
squares values of the parameters will be recorded by 
SEAC (along with the deviations between theory 
and experiment). This augmentation reduces the 
capacity of the code when used with the high-speed 
memory; with 18 adjustable parameters, a maximum 
of 22 observations could be fitted, whereas with 9 
adjustable parameters, 63 observations could be 
fitted. However, an alternative slightly slower code 
is available if external tapes are used, and in this the 
only limitation is that the sum of the number of 
observations and the number of parameters be less 
than 256. 

The Ta n problem was worked in an approxima- 
tion where first 8 and then 9 parameters were ad- 
justed (i. e., results were obtained with and without 
the final parameter a included). There were 39 ob- 
served values, but the number of these was reduced 
in successive calculations. As a check, both the 
observed values and the original deviations from 
these values were fitted. It required about 20 min 
of machine time for each of these calculations as com- 
pared to about 9 hr required in an earlier hand com- 
putation, where the normal equations were first 
obtained and then solved by elimination. Thp 
suits are described in section 7.3. 

" P. Davis and P. Rabinowitz, J. Assoc. Comp. Mach. 1, 183 (1954). 
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5. Data Preparation — Matrix Inputs 

A code has been set up that will generate all 
matrices with a single input of the parameter values. 
It uses the internal memory and can generate 
matrices of order up to 38X38. Checks in the code 
insure that the matrix generated is correct on wire 
and ready as an input for the coding described in 
section 3. For the Ta n work, all matrices could 
be generated in 30 to 60 min of machine time. 

When a problem is first started, it is necessary to 
k type many coefficient matrices (20 for the Ta n 
work) for each final matrix generated, and it would 
seem simpler to evaluate the elements (1) with a 
desk computer and type the final matrix directly. 15 
However, only the nonzero coefficients need be 
typed, and as there are usually very few of these, 
the total work involved is about the same. 16 Once 
the coefficient matrices are typed, final matrices can 
be generated for as many sets of parameter values as 
desired. In any one problem, matrices must be set 
up with at least two sets of parameter values, so an 
appreciable amount of time will be saved. 17 The 
same coefficient matrices are used to generate the 
linear formulas (2) with the coding described in 
section 6, and in relation to this work, the time 
required to type the coefficient matrices is a negligible 
factor. 

The coefficient matrices used with this code must 
have no more than 56 nonzero elements, but larger 
coefficient matrices can be inserted piecemeal by 
associating extra parameters (which will be given 
the same value) with each part. To allow for this 
possibility, the code has been written for 24 param- 
eters, which is probably a larger number than would 
ever be needed. The largest coefficient matrix 
encountered in the Ta n problem had only 36 
nonzero elements. 

6. Data Preparation — Least-Squares Inputs 

A code has been set up that will generate the co- 
efficients, formula (3), that enter into the linear 
formula for the eigenvalue, formula (2). Instruc- 
tions that determine the number of independent 
parameters desired in the linear formula outputs are 
first fed in; frequently, as here, parameters having 
the same symbol are regarded as identical in all 
configurations in which they occur, but other 
approximations are also used. All vectors for which 
linear formulas are desired are then fed into the 
internal memory ; there is space for 582 vector com- 

1 5 The time required to set up the code was naturally far in excess of the time 
needed to do the work with hand computers. The expenditure of this time is 
justified only because it is expected that the same coding will be used subse- 
quently in many similar problems. 

16 The last four digits of each coefficient are used to specify the position in the 
matrix; this leaves seven digits for the coefficient itself. In the present work 
only four significant figures were retained after the decimal point, because a 
result theoretically correct to within 1 kayser was considered adequate with a 
mean deviation between theory and experiment more than a hundred times as 
great. 

17 The same set of coefficient matrices can often be used to generate final matrices 
for other spectra that have the same, or related, configurations. The coefficient 
matrices used in the Ta n work also apply to the configurations d 8 , d 7 s, and rf 6 s 2 , 
and they have been used to generate the matrices for the first spectrum of ruthe- 
nium (Ru i) as a start toward assisting with the analysis of this spectrum. 



ponents (all vectors of a 24X24 matrix). 18 The 
coefficient matrices are then fed in, and for each 
one fed in (or for each group of consecutive matrices, 
depending on the instructions), as many coefficients 
of linear formulas print out as there are vectors in 
the vector input (the maximum number of possible 
outputs is 24). For the Ta n work, the 39 linear 
formulas in nine parameters were computed in less 
than 60 min of machine time; originally this work 
was carried out with a desk computer, and the time 
required was of the order of 60 hr. 

The coding just described groups coefficients of 
the same parameter only if they are in the same 
matrix. For use in the least-squares code of section 
4, coefficients of the same parameter must be grouped 
together without regard to the J-value of the matrix 
in which they occur. The second part of this code 
sorts out the preceding coefficient outputs. A 
remaining uncorrected deficiency is that it does not 
give an output ready to be used as an input for the 
orthonormalization code of section 4; the augment- 
ing and check sums were inadvertently overlooked, 
so the results had to be outscribed, amplified, and 
reinscribed to get the correct input. This sorting 
can be carried out only for problems where the 
number of coefficients, formula (3), is less than 672. 
With 15 parameters (which may be the least number 
that can give good agreement in Ta n) only 44 
observations can be handled (and 39 of 81 possible 
levels have already been observed in Ta n). This 
limitation is a very difficult one to overcome without 
altering the construction of the computer itself. 
However, it is always possible to combine two 
partial sortings and set up by hand the final input 
for the least-squares code. 

It would be convenient at this point to carry out a 
few manipulations with the sorted linear formulas, 
but the potentialities for doing this have not yet 
been fully explored. However, supplemental coding 
has been made that will evaluate the linear formulas, 
and also the differences between the formulas and a 
set of observed values, for as many sets of parameters 
as desired. 

7. Results for Ta II 

7. 1 . Preliminary Estimates of Parameters 

The spectrum of Ta n is very complicated because 
of the presence of the three overlapping configura- 
tions of 5rf 4 , 5d? 6s, and 5d 2 6s 2 . The analysis of the 
experimental data 19 is far from complete, as only 
39 of 81 possible levels in these configurations have 
been identified. «7-values are assigned to all identified 
levels, but the configuration and XS-coupling desig- 

is With some amplification, the coding could be used to obtain nondiagonal 
elements of the transformed coefficient matrices (as well as the diagonal elements, 
formula (3)); it is estimated that matrices of order up to 18 X 18 might be trans- 
formed with this amplified code. Such a code would be useful in solving the 
problem referred to in footnote 6 because it is sometimes easiest to set up coeffi- 
cient matrices in a scheme that does not correspond to observed behavior, and 
then to use a simple transformation to arrive at a labeling that is more nearly 
approached in nature. 

i* C. C. Kiess, G. R. Harrison, and W.J.Hitchcock, J. Research NBS 44. 245 
(1950) RP2075. 
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nations have been made for only 26 of them. How- 
ever, the configuration assignments of six of the 
lowest levels have been uncertain. Kiess (see foot- 
note 19) had assigned these levels to the 3 F and 3 P 
terms of the 5d 2 6s 2 configuration, but this assign- 
ment put the 3 P term below the 3 F term, whereas 
it would be expected to be much higher on the basis 
of a simple iS-coupling theory. It would have been 
reasonable to assign the 3 F term to the bd? 6s con- 
figuration, but this would leave the 3 F term in 
3d 2 6s 2 still unidentified and force one to conclude 
that the lowest levels of this spectrum, corresponding 
to the bd 2 6s 2 3 F term, were completely unknown. 20 In 
any case, it was certain that configuration and LS- 
coupling assignments would be approximations at 
best, so that it would be very hard to decide what 
experimental values should be correlated with 
theory in estimating preliminary values for the 
parameters. 

Because of this, the theory was first applied to the 
experimental values 21 in the simpler, almost com- 
pletely analyzed spectrum of singly-ionized hafnium 
(Hf n), to get some idea of the values that would be 
expected for the parameters in spectra with 5d- 
electrons; no calculations had been published for 
these spectra at that time. The Hf n parameters 
were not evaluated by least-squares, but the values 
obtained agreed well with those recently published 
by Gehatiah (see footnote 3). 22 

Preliminary parameter values were then estimated 
for Ta ii by adjusting the Hf n values to get rough 
agreement with a few of the best identified levels of 
the Ta n spectrum. The magnitude of the read- 
justment was kept reasonably consistent with what 
would be expected from the known behavior of the 
parameters in spectra with 3d- and 4rf-electrons. 
However, some of the parameters were poorly ad- 
justed, and acceptable preliminary values of the 



2° In commenting on this, Kiess noted that too few strong transitions were left 
unidentified in his line list to give much weight to this possibility (private 
communication). The assignment made by Kiess was originally made also by 
P. F. A. Klinkenberg, G. J. van den Berg, and J. C. van den Bosch, Physica 16, 
861 (1950), but they later favored the second possibility that the lowest levels had 
not yet been identified (Physica 17, 167 (1951)); theory was used to demonstrate 
the validity of the second alternative by van den Berg (Ph. D. Thesis, Uni- 
versity of Amsterdam, 1951 (unpublished)) . His calculation is much too simpli- 
fied to decide the point, however, and the calculation given here shows conclu- 
sively that the first alternative is the correct one. 

2i W. F. Meggers and B. F. Scribner, J. Research NBS 13, 625 (1934) RP732. 

22 Gehatiah's parameter values for Hf n are: 

£=435, C= 1,530, G 2 = 2,818, H 2 = 482, and r = 1,336. 

Analogy with results for 3d- and 4d-electrons suggests that the Ta n parameters 
should more nearly equal to the parameters of Ta i or W i (excepting the spin-orbit 
parameter f for the latter spectrum) than to those of Hf n. Some recent unpub- 
lished evaluations of parameters in these two first spectra are interesting. Racah 
has given the final parameter values that he obtained in the 5rf 4 6s 2 and 5d 5 65 con- 
figurations of W 1: 

5=371, C= 1,900, ^2=2,700, H 2 =434, r = 2,089 and a =46. 

His calculations are briefly described in footnote 7. The Ta 1 calculations are 
for the 5^3 6s 2 configuration and neglect configuration interaction; they have 
been carried out with and without the L (Z-f-1) correction (see R. E. Trees, 
Phys. Rev. 92, 308 (1953) and also T. Kamei, Phys. Rev. 99, 789 (1955): 

JB=345, C = 1,289, f = 1,657, a = 112. 

B = 278, C = 2,100, f = 1,650. 

The neglect of configuration interaction in the Ta 1 calculations is probably not 
fully justified, which may account for the value of a. being larger than expected. 



parameters were still obtained, so it is likely that 
this readjustment is not critical. Some details of 
the procedure will be given, however, mainly to 
bring out significant features of the parameter 
behavior. The preliminary parameter values are 
given in calculation I of table 1, and the final least- 
squares-adjusted values of the parameters are given 
in calculation II. 

The preliminary values for parameters B and O 
were assumed to be 20 percent greater than the 
values in Hf n, although an increase of only about 15 
percent would be expected from analogy with 3d- 
electron spectra. Calculation II shows not only 
that this overestimated B and C, but also what is 
more important, that the final ratio of B to C in 
Ta ii is about 8 percent smaller than the ratio in 
Hf ii (it is 25 percent greater than the ratio for the 
W i parameters (see footnote 22). This ratio has a 
more constant value in spectra with 3d- and 4d- 
electrons (see footnote 22). The ratio may vary in 
spectra with 5rf-electrons 3 but it is also possible that 
a nearly constant ratio would be obtained in these 
spectra if a higher-order approximation were used, 
as suggested at the end of this paper. 

The preliminary estimate 6? 2 = 2,000 was made to 
favor the assumption that the lowest identified 3 F 
term belonged to the bd? 6s configuration, even 
though this assumption was not considered correct. 
Calculation II shows that G 2 should have been given 
the value found in Hf ii. Analogy with 3 (/-electron 
spectra indicates that G 2 should be slightly larger in 
Ta ii than in Hf n, but the increase is much less than 
expected for B and O. 

The value f=2,100 is a considerable overestimate; 
calculation II indicates that a value slightly larger 
than the value applicable in Ta i should have been 
used (see footnote 22), which is what would be ex- 
pected from analogy with 3d-electron spectra. 

To estimate the interaction parameter H 2 , the 3 F 
and 3 P matrices were set up in L$-coupling by using 
all the preceding preliminary parameter evaluations. 
The lowest terms of these matrices were expected to 
be very sensitive to the value of the parameter H 2 , 
so that errors in estimating the other parameters 
would not be too important. A value of 500 was 
assumed for H 2 , as it was expected that the parameter 
would have a value in Ta n that would nearly equal 
that in Hf n. Analogy with the behavior shown in 
the 3d-electron spectra would lead one to expect that 
the value in Ta n might be smaller than in Hf n, but 
as this integral depends strongly on the overlapping 
of d- and s-wave functions, the analogy would not be 
very reliable. The eigenvalues of the two 5X5 
matrices were obtained on SEAC with the coding de- 
scribed in section 3. The lowest 3 P term was found 
to be 800 K above the 3 F term, instead of 1800 K 
below as observed. It was decided that an extra 10- 
percent increase in the value of H 2 would be justified, 
and that this, together with the inclusion of second- 
order spin-orbit effects, would explain the observa- 
tions with reasonable parameter values. 
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Table 1. Configurations d 2 s 2 , d 3 s, and d* of Ta n in intermediate coupling 



Levels are identified by their ./-values and by the label number, which is the order in energy sequence. 

Calculation I: Solutions of matrices for preliminary values of the parameters. 

Calculation II: Solution from linear formulas; parameters evaluated by least-squares, with the levels 18500 and 23295 omitted. 













Calculation I 




Calculation II 


J-value 


Label 


Observed 


Obs. g 
























energy 






















Calc. energy 


Calc— obs. 


Calc. g 


Calc. energy 


Calc— obs. 





1 
2 
3 
4 
5 
6 

( 1 


4125 
12601 




3798 
11591 
15240 
23317 
25860 
33923 

-292 


-327 
-1010 




4046 
12547 


-79 
-54 










23381 




-64 




23084 


-297 





















0.000 


-292 


0.058 


-53 


-.53 




2 


5331 


1.550 


5316 


-15 


1.537 


5216 


-115 




3 


10713 


2.374 


11351 


638 


2.128 


10982 


269 




4 


13475 


1.510 


13111 


-364 


1.654 


13331 


-144 


1 


5 
6 


14628 
17375 


0.850 
1.171 


13651 
17014 


-977 
-361 


1.119 
1.007 


14156 
17175 


-472 
-200 






8 
I 9 

/ 1 






22412 
26426 
30446 

853 


































1031 


1.008 


-178 


1.009 


1139 


108 




2 


3180 


0.750 


2607 


-573 


0. 767 


3532 


352 




3 


5658 


1.340 


5045 


-613 


1.126 


5872 


214 




4 


9690 


1. 063 


7865 


-1825 


1.039 


9469 


-221 




5 


11875 


1.48 


11706 


-169 


1.210 


11653 


-222 




6 


13560 


1.11 


13892 


332 


1.495 


13246 


-314 


2 


7 


14494 


1.472 


14794 


300 


1.485 


14375 


-120 




8 


16424 


1.405 


16576 


152 


1.392 


17227 


803 




9 


17168 


1.211 


17966 


798 


1.289 


17886 


718 




10 


18500 


1.462 


22389 


3889 


1.094 


22836 


(4335) 




11 


22928 


0.700 


23128 


200 


0.732 


22866 


-63 




12 
\ 13 

/ 1 






28091 
30795 

2980 






















2642 


1.250 


338 


1.244 


2663 


21 




2 


6831 


1.098 


6598 


-233 


1. 080 


0K20 


-11 




3 


11767 


0.915 


9876 


-1891 


0. 996 


11790 


23 




4 


12435 


1.614 


13679 


1244 


1.467 


12730 


295 




5 


14581 


1.004 


13322 


-1259 


1.048 


13293 


-1288 


3 


? 


15726 


1. 476 


16197 


471 


1.460 


I5S9S 


172 




23295 


1.120 


19144 


-4151 


1.348 


18363 


(-4932) 




8 


23620 


1. 076 


23976 


356 


1. 080 


23464 


-156 




9 


24870 


0. 995 


24563 


-307 


1. 006 


24666 


-214 




10 


26829 


0.855 


27699 


870 


0.852 


27054 


225 




V 11 
/ 1 






32217 
5068 










4416 


1. 350 


652 


1.329 


4434 


18 




2 


9746 


1.225 


9421 


-325 


1.174 


9651 


-95 




3 
4 






11206 
12375 




1.110 
1.061 






12705 


1.021 


-330 


12685 
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7.2. Preliminary Calculation of Energy Levels 

By using these preliminary values for the para- 
meters, matrices were set up for the Tan problem, as 
outlined in section 2. (The coding in section 5 was 
not then available.) The eigenvalues and eigen- 
vectors of these matrices were then calculated on 
SEAC, as described in section 3. The results are 



given in calculation I of table 1 . Partly because the 
eigenvectors obtained in this calculation are based on 
preliminary values for the parameters, no attempt 
has been made to give the conventional Z#-coupling 
designation for the levels, and the levels have been 
identified by the J- values alone; the assignments 
given to some of these levels by Kiess (see footnote 
19) are, however, roughly confirmed. 
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Although the preliminary parameters were esti- 
mated without making any effort to assign the 
majority of the observed levels to any theoretical 
state, a good correlation can be established between 
observed energy levels and calculated eigenvalues, 
and between the corresponding observed and cal- 
culated ^-values. There are, however, three observed 
levels with J=2 (16424, 17168, and 18500) in a 
region where there are only two calculated (16576, 
17966) ; the (/-values indicate that the level which can- 
not be correlated with theory is the one observed at 
18500. Similarly, there are three observed levels 
with J=3 (23295, 23620, and 24870) in a region 
where two are calculated (23976 and 24563) ; in this 
case the ^-values indicate that the level which cannot 
be correlated with theory is the one observed at 
23295. The levels with J=3 observed at 12435 and 
14581 were correlated with the respective calculated 
values 13679 and 13322. This is the only instance 
where the calculated and observed values were not 
correlated in the same sequence with respect to 
energy. This correlation was made when it was 
found with better parameter estimates that the linear 
formulas for the eigenvalues indicated that their 
energies would cross. This correlation could have 
been made more easily by noting that it makes the 
^r-values agree much better, but this criterion would 
also call for reversing the correlations of the fifth and 
sixth levels with J=2, and this was not done. With- 
in the limits of accuracy of this preliminary calcula- 
tion, the only other correlations that seem uncertain 
are the ones connecting the level observed at 12705 
with the level calculated at 12375, rather than the 
ones calculated at 11206 or 14675, and also the cor- 
relation of 17231 with 17597 rather than with 18354. 
Both these correlations worked well in subsequent 
least-squares calculations. However, the least- 
squares calculations show that agreement would be 
improved if the level observed at 16424 were disre- 
garded, rather than the one at 18500. It is likely that 
other possibilities still remain, that might give better 
agreement. 

7.3. Least-Squares Estimates of Parameters 

Linear formulas for the eigenvalues were set up by 
hand (the coding of section 6 was not then available) 
in an approximation with nine adjustable parameters ; 
the parameter a, which was neglected for the pre- 
liminary estimates, is the additional parameter. 
Least-squares calculations were then carried out on 
SEAC, with and without the parameter a, as de- 
scribed in section 4. These calculations will be 
described briefly; the parameter values that were 
obtained are given in successive columns of table 2, 
labeled to correspond to the descriptions. 

(I) The parameters were determined by using all 
39 observed levels with the correspondence between 
theory and observation the same as that already 
given in calculation I of table 1. The purpose was to 
verify, what was already fairly certain, that no 
readjustment of the parameters would give a good fit 
for all 39 observations. In this calculation the level 



Table 2. Least-squares evaluations of the parameters of 
d 2 s 2 , d 3 s, and d 4 of Tan 

The parameters A determine the centers of gravity of the three configurations; 
B, C, f, and a would be specified independently in all three configurations in 
higher approximations than the one given here. The interaction parameter, 
G2(d i —d 2 s 2 ), has been taken equal to the exchange integral G2(d 3 s) and desig- 
nated Gr, the two interaction integrals //2(d 4 -d 3 s) and H2(d* s—d^s 2 ) are taken 
equal and designated H2. 

The first of a pair of determinations grouped together omits the L(L-\-\) cor- 
rection, the second includes it. A correspondingly numbered section of the text 
specifies the observations included in making each pair of calculation. 
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1947 
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539 
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26830 
479 
1824 
2924 
1843 
521 


20388 
13957 
26356 
447 
2012 
2918 
1874 
544 
-35 
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at 18500 had a deviation of 3817, and the deviation 
for the level at 23295 was -4294. The level at 
12831 had the next largest deviation of 1316, which 
is an order of magnitude smaller. The mean devi- 
ation was ±1048; the mean deviation, exclusive of 
the two badly fitting levels, was ±515. 

The inclusion of the L(L-\-l) correction produced 
a negligible improvement in the mean deviations, 
and the correction was negative (a =—35), whereas 
all the evidence so far accumulated indicates that it 
should be positive. With the Z(Z + 1) correction the 
mean deviation for all levels was reduced to ±1,038, 
and the deviation, excluding the two badly fitting 
levels, was ±510. 

(II) The levels at 18500 and 23295 were omitted 
from the least-squares calculation; the deviations for 
this calculation are given in calculation II of table 1. 
The mean deviation of this calculation is ±371: As 
this is considerably less than the value ±515 for 
calculation (I), described in the previous paragraph, 
it indicates that the parameters of the latter calcu- 
lation were influenced by the inclusion of the two 
badly fitting levels. Comparison of the parameters 
themselves also indicates this. In this calculation 
there are three deviations (for the levels 16424, 
14581, and 12831) that exceed twice the mean devi- 
ation of all levels fitted. 

Again, when the L(L+1) correction is included, 
the mean deviation is reduced only slightly (to the 
value ±364) and the correction is negative (a= — 19). 
Most of the improvement produced by including the 
Z(Z + 1) correction can be attributed to the slightly 
better agreement obtained for the badly fitting level 
12831, which has its deviation reduced to 858; if the 
mean deviation is computed without this level the 
value ±371 would be replaced by ±341, and the 
value ±364 replaced by ±340. The inclusion of 
the L(L-\-l) correction is therefore not justified. 

As calculations III and IV indicate, the inclusion 
of three badly fitting levels has not had much effect 
on the agreement obtained for the other levels. 
Also, there seems to be no possibility of including 
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the L(L-\-\) correction and getting an appreciable 
improvement. The best calculation that can be 
offered at this time is that in calculation II of table 1. 
Possibilities for improving this calculation will be 
given in the discussion. 

(III) The levels L8500, 23295, and 12831 were 
omitted from the least-squares calculation. The 
mean deviation was ±335, which is only a little less 
than the value ±341 obtained in (II). Inclusion 
of the L(L-\-l) correction produced negligible im- 
provement, as the deviation was reduced only to 
±333; the correction was still negative (a= — 10). 
It would have been difficult to justify the omission 
of the level 12831 if this calculation had given much 
better agreement. 

(IV) Eleven levels were omitted in this calculation. 
These omissions consist of the two pairs, 18500 and 
22928, and 12435 and 14581, omitted partly because 
the calculated eigenvalues were close and it seemed 
possible, therefore, that the linear formulas might 
not apply; the two triplet groups 23295, 23620, and 
24870, and 16424, 17168, and 18500, which were 
omitted because of uncertainties in the correspond- 
ence between theory and experiment; and the levels 
23381 and 12831, omitted because they had rather 
large deviations that tended to be improved by a 
negative L(L-\-l) correction. The mean deviation 
of this calculation was ±170 (the mean deviation of 
the same terms in calculation (II) is ±200). The 
only large deviation was —491 for the level 14628. 
With the i(Z+l) correction included, the mean 
deviation was reduced to ±151, due mainly to a 
decrease in the deviation for the level 14628 to 
— 306. Again the sign of the correction was nega- 
tive (a =-35). 

8. Discussion 

Results obtained in spectra with 3d- and 4r/-elec- 
trons indicate that when the L(L-\-l) correction is 
used the mean deviation between theory and experi- 
ment will generally be less than 200 K, and that the 
sign of the correction will be positive. 23 Calculation 
(II) of table 1 should therefore be improved because 

23 Racah (see footnote 7) reports a mean error of 411 without the L(L+1) cor- 
rection, and an error of 229 with the L(L+l) correction in W i; as indicated in 
footnote 22, the sign of the L(L-\-l) correction is positive. By coincidence, the 
mean error of calculation II, table 1 is 412 K, nearly identical with the corre- 
sponding W i value; however, no appreciable improvement was obtained by the 
use of the L(L+\). correction In Hf n, Gehatiah (see footnote 3) obtains a mean 
error of 257 without the use of this correction. 



a better significant agreement between theory and 

experiment can probably be obtained. 

It is fortunate that the experimental analysis of 
the Ta n spectrum will be extended by members of 
the NBS Spectroscopy Section. The calculation 
given is already accurate enough to furnish con- 
siderable help. This analysis may, in turn, help 
solve the problem of the two excess levels pointed 
out previously. In a preliminary survey, the switch- 
ing of J-values has been ruled out, as none of the 
levels concerned are poorly identified. The identi- 
fication of new levels will help to eliminate some of 
the arbitrariness involved in correlating the observed 
and calculated values. With more observations a 
calculation with extra adjustable parameters (sug- 
gested below) will also be better justified. 

Eigenvalues and eigenvectors of the matrices set 
up with least-squares values of the parameters will 
be obtained on SEAC after the coding now under- 
way is completed (see footnote 13). The g- values 
obtained should be much more accurate than those 
already calculated, and if this is so, the uncertainties 
in making correlations may be resolved. Knowl- 
edge of the exact eigenvalues will show whether or 
not any of the large errors of the calculation are due 
to limitations on the validity of linear formulas for 
the eigenvalues. The linear formulas obtained from 
these new solutions should be accurate enough so 
that eigenvalues can be redetermined from them for 
any additional small changes that are made in the 
parameters, without setting up new matrices. 

The agreement between theory and experiment 
may also be improved by using more adjustable 
parameters; B, C, and f can be adjusted independ- 
ently in all three configurations, thus increasing the 
number of adjustable parameters from 9 to 15. 
Theoretically, this should be done, as the effect of 
differences in these parameters is likely to have the 
same order of magnitude as the L{L-\-\) correction. 24 
The difficulty of carrying out least-squares calcula- 
tions with this number of parameters is great, even 
with SEAC. Some coding improvements and more 
experience in the behavior of the codes with this 
type of problem are required. 



24 See section 4.2 of the reference in footnote 2, 
footnote 7. 



or section 2 of the reference in 



Washington, September 9, 1955. 
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